home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
CD ROM Paradise Collection 4
/
CD ROM Paradise Collection 4 1995 Nov.iso
/
os2
/
adaptor.zip
/
ADAPT.ZIP
/
adaptor
/
examples
/
purdue
/
prob01.fcm
next >
Wrap
Text File
|
1993-06-26
|
2KB
|
72 lines
PROGRAM PROB01
C
C PROBLEM 1
C
C REFERENCE: PROBLEMS TO TEST PARALLEL AND VECTOR LANGUAGES
C CSD-TR 516, COMPUTER SCIENCE, PURDUE UNIVERSITY
C JOHN R. RICE, MAY 1, 1985
C
C REVISED BY JOHN R. RICE AND J. JING, OCT. 1, 1990
C
C
C *************************************************
C * Adapted for FORTRAN D benchmarking *
C * by T. HAUPT (haupt@sccs.npac.syr.edu) *
C * *
C * Northeast Parallel Architectures Center *
C * at Syracuse University, Syracuse, NY, USA *
C *************************************************
C
C
C VERSION SIMD/CM2-1.00
C ==================================================
C
INCLUDE '/usr/include/cm/paris-configuration-fort.h'
INTEGER KASES,K,NK
PARAMETER (KASES=6)
INTEGER N(KASES)
CMF$ LAYOUT N (:serial)
REAL TN
REAL A,B,H
DATA A,B / 0.0,1.0 /
DATA N / 8192, 65536,
X 131072, 262144, 524288, 1048576/
C
C LOOP OVER KASES
C
DO K = 1, KASES
NK=N(K)
CALL CM_TIMER_CLEAR(0)
CALL CM_TIMER_START(0)
DO MANY=1,100
H = (B-A)/(NK-1)
CALL DOIT(NK-1,A,H,TN)
TN = H*((EXP(A)+EXP(B))/2.+TN)
ENDDO
CALL CM_TIMER_STOP(0)
PRINT 30,A,B,NK
30 FORMAT ('PROBLEM 1 WITH A,B,N = ',2F10.3,5X,I8)
PRINT 40,TN
40 FORMAT ('GIVES TN =',F10.6)
CALL CM_TIMER_PRINT(0)
ENDDO
END
SUBROUTINE DOIT(N,A,H,FSUM)
INTEGER N
REAL A,H,FSUM
REAL R(N)
R = EXP (A + H * [1:N])
FSUM=SUM(R)
END